Proximity Effect and Multiple Andreev Reflections in Chaotic Josephson junctions 
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We study the dc-current transport in a voltage biased 
superconductor-chaotic dot-superconductor junction with an 
induced proximity effect(PE) in the dot. It is found that for 
a Thouless energy Eth of the dot smaller than the supercon- 
ducting energy gap A, the PE is manifested as peaks in the 
differential conductance at voltages of order Eth away from 
the even subharmonic gap structures eV w 2(A ± Eth) /2n. 
These peaks are insensitive to temperatures kT <C A but are 
suppressed by a weak magnetic field. The current for sup- 
pressed PE is independent of Eth and magnetic field and is 
shown to be given by the Octavio-Tinkham-Blonder-Klapwijk 
theory.™ 

Over the last decades, there has been a large interest in 
various aspects of the proximity effect (PE) in mesoscopic 
normal conductor-superconductor (NS) systems. The PE 
can on a microscopic level be viewed as correlations on 
the scale of the Thouless energy Eth between electrons 
and holes in the normal conductor, induced via Andreev 
reflections at the NS-interface. Well known manifesta- 
tions of the PE are the bias conductance anomalycl in 
diffusive NS-junctions and the induced gap in the spec- 
trum of diffusive^ or chaoticQ SNS- and NS-junctions. 

There is however no complete theory of the PE in volt- 
age biased SNS-j unctions. The reason is that electrons 
and holes in the .Darmal conductor undergo multiple An- 
dreev reflectionarEil (MAR), giving rise to complicated 
correlation effects and staong nonequilibrium in the nor- 
mal part of the junction.tj The fingerprint of MAR trans- 
port is subharmonic gap structures (SGS) in the current- 
voltage characteristics at eV = 2A/n. 

So far only some limiting cases have been considered. 
In short junctions, where Etk is much larger than 
superconducting energy gap A, a coherent MAR-theor 
has been developed which fully incorporates the PE. The 
theory describes to^arge accuracy experiments in atomic 
and has also been applied toj-short 
and disordered tunnel barriers li!J For 
junctions in the short limit, the PE modifies the shape 
but not the position of the SGS, as well as gives rise 
to ac- Josephson effect. Other cases studied are long, 
Etk *C A, diffusive SNS-junction with no proximity cou- 
pling between the superconductorsO and SNS-j unctions 
consisting of t\sia weakly coupled proximity NS-junctions 
in equilibrium .113 

Taken together, present theories do not give an answer 
to the general question of the role of the PE in SNS- 
junctions where the Thouless energy is comparable to the 
superconducting energy gap, Etk ~ A. Moreover, recent 
experiments on SNS-junctions in this regimetJ show qual- 



DrJ 



size point contact 
diffusive junctions 



itatively new features: splitting of the SGS conductance 
peaks at even subharmonics eV = A and A/2. 

In this Letter we study current transport in 
a superconductor-chaotic dot-superconductor (S-dot-Si 
junction. By employing the scattering theory of MAREJ 
and using random matrix theoryO to describe the statis- 
tical properties of the energy- and magnetic field depen- 
dent scattering matrix of the dot, we are able to calculate 
the dc-current for arbitrary ratio between Eth and A as 
well as for arbitrary strength ofj-the PE due to variation 
of the magnetic field in the dot.Q 

The main result of the paper is the explanation of the 
splitting of the even SGS conductance peaks, i.e the ad- 
ditional peaks at eV rj 2(A ± ETh)/2n, in terms of an 
induced PE in the dot. The PE conductance peaks are 
suppressed by a weak magnetic field but are insensitive 
to temperatures kT <C A. 
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FIG. 1. a) Schematic picture of the junction. The shaded 
region is the chaotic dot and the black bars denote tunnel bar- 
riers with transparency Fl and Tr. b) A MAR-trajectory in 
energy-position space for quasiparticles injected at energy E 
from left. Filled (empty) chaotic dots denote electrons (holes) 
scattering. Scattering processes in which quasiparticles are 
Andreev reflected at energies Eth away from the Fermi en- 
ergy have an enhanced amplitude, giving rise to additional 
peaks in the conductance at eV w 2(A ± i?Th)/2n. 

The junction is shown schematically in Fig.|f|: A two- 
dimensional quantum dot is coupled to two supercon- 
ducting electrodes via quantum point contacts support- 
ing N transverse modes each. The contacts contain tun- 
nel barriers with mode-independent transmission proba- 
bilities Tl, Tr ^> 1/N. It is assumed that the quasiparti- 
cle dwell time in the dot, U/Etk: is much smaller than the 
inelastic scattering time; here Etk = N(Yl + Fr)S/(2it) 
and 6 is mean level spacing of the dot. In this case the 
transport through the junction can be characterized by 
the Andreev reflection amplitudes at the normal lead- 
insulator-superconductor (NIS) interfaces and the scat- 
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tering matrix of the dot. We consider the case where the 
classical motion in the dot is chaotic on time scales longer 
than the ergodic time r erg ~ L/vf {I 3> L) or L 2 /(Ivf) 
(I -C L), where L and I are the linear dimension and the 
mean free path of the dot. The ergodic time is assumed to 
be smaller than the quasiparticle dwell time and the in- 
verse superconducting gap_^r erg <C Ti/Erh, h/A. In this 
case random matrix theoryE3 can be used to describe the 
scattering properties of the dot. 

The scattering matrix S can be written in terms of the 
Hamiltonian H (dimension M) of the closed dot as 




= 1 - 2mW\E -H + iirWW^^W, 



(1) 

where r, t, r' and t' are the N x N reflection and trans- 
mission matrices and W describes the coupling of the dot 
to the leads, with W nm = 8 nm (MS) 1 / 2 /Tt. The Hamil- 
tonian H is described by a random Hermitian matrix 
H = Hq + ijHi, where H and Hi are real symmet- 
ric and anti-symmetric matrices respectively, indepen- 
dently distributed with the same Gaussian distribution, 
P(H 0{1) ) cxcxp[-7r 2 (l+ 7 2 )tr( J ff 0(1) iJ T (1) )/(4M ( 5 2 )]. The 
parameter 7 is related to the magnetic flux in the dot as 
$ ~ 7$o(-^T" erff /7i) 1 / 2 ,where $0 = h/e is the flux quan- 
tum. A magnetic flux <& c ~ &o{ T ergETh / h) 1 ^ 2 effectively 
breaks time reversal symmetry in the dot. 

The current is calculated within a scattering approach 
for the Bogoliubov-de Gennes equationstLa and is ex- 
pressed in terms of the scattering state wavefunctions 
as 1= (e/h) J dy^hxi^lid^a/dxDfo, where a = 
{e/h, E, L/R,j} labels the scattering state (electron/hole 
like quasiparticle injected at energy E from the left/right 
superconducting reservoir, in transverse mode j) and fo 
is the equilibrium Fermi distribution of the reservoirs. 
The energies are measured relative to the Fermi energy 
in the superconductors and the magnetic field in the con- 
tacts is assumed to be zero. 

The acceleration of the injected quasiparticles due to 
the applied voltage V gives a scattering state wavefunc- 
tion which is a superposition of electron and hole states 
at different energieaij E„ = E+2neV. For quasiparticles 
incident from the left superconductor, the wavefunction 
on the normal side of the left NIS-interface has the 

N 

*L = E $ ™(y)E e " 4iS2 " t/ri 

m— 1 n 
' ( r e '+ p ik l,,2n x + r e - p- ik m,2n x \/(he U/2 

(J*,+ P ik^ 2n x, h - p ~ik^ 2n x\j (k h \l/2 

where m is transverse mode index. At the right i 
face the electron/hole energies are shifted by ±eV and 
the wave function is given by multiplying $i by 
exp(— ia z eVt/h) and substituting 2n — > 2n + 1 (a z is 
the Pauli matrix in electron-hole space). The wave vec- 
tor km h n — [{2m/h 2 )(EF — e m iEn)) 1 / 2 , where e m is the 



transverse mode energy and +/— corresponds to elec- 
trons/holes. The vector potential enters only the trans- 
verse wavefunctions & m , normalized for each mode to 
carry the same current. The scattering at the dot con- 
nects the electron and hole wavefunction coefficients as 

(£)-*(£)•(£)-*(&)• <" 

where we have introduced the vector notation c^' + = 
[c?',;t. and Si = S(E n ) and S~ = S*(-E n ). At 

the left NIS-interface, the scattering is described by An- 
dreev and normal reflection and transmission amplitudes 
for electrons (e) and holes (h), af/ h — a e / h {E n ) and sim- 
ilarly bn , Cn h and dt/ h , given in Ref . |f6|. This gives the 
connection between wavefunction coefficients as 

(S) -(5 $)(£)♦**»(*)• (4) 

where the source term (oc S n oS m j) describes electron 
quasiparticle injection. The coefficients at the right NIS- 
interface are connected in a similar way. The other scat- 
tering states are constructed analogously. 

In the short dwell time regime, Exh A, the scat- 
tering matrix is independent of-energy on the scale of A, 
and the current can be writtenBli3 as a sum of the single 
mode currentso, with different transmission eigenvalues 
D m (eigenvalues of the matrix product tv). The en- 
semble averaged current (I) is then found via an integral-, 
tion over transmission eigenvalues with the distributiontll 
pp) = N/n[T(2 - T)]/([r 2 + 4D(1 - Y)\JD{\-D)). 
The current voltage characteristics if*-|Fig. ga show SGS 
at eV — 2A/n and an excess currentta for all T [the nor- 
mal state conductance G N = {2e 2 /h)NT L T R /{T L +T R )]. 




For a longer quasiparticle dwell time, Etk % A, the 
ensemble averaged current is calculated numerically by 
generating a large number of HamiltoniansO The cur- 
rent voltage characteristics is shown in Fig. 0b for 
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FIG. 2. The current- voltage characteristics for a symmetric 

junction Y = Tl — with a temperature kT < A. a) In 
the short dwell time regime, Eth 2> A, the current show SGS 
at eV — 2A/n. b) In the intermediate dwell time regime, 
Eth = 0.4A, the current is shown for $ = (solid) and 
$ 3> $ c (dashed) and N = 10. The SGS is less pronounced 
and the current for low transparencies is smaller compared to 
the short dwell time regime. 
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Eth = 0.4A. Compared to the short dwell time regime 
in Fig ||a, the SGS is less pronounced and the current for 
low barrier transparencies is reduced for $ = and even 
further for $ > $ c . 

The current voltage characteristics can be studied in 
detail by considering the conductance, dl/dV, shown 
in Fig. H (parameters chosen to clearly display PE- 
features). The conductance with suppressed PE, $ 3> 
$ c , shows SGS at eV = 2A/n. The PE is manifested 
as an enhancement of the SGS at eV — A but also as 
additional peaksnJ in the conductance on both sides of 
eV = A. It is found by varying Exh that the positions of 
the additional peaks are given by eV « A ± Exh ■ More- 
over, the position and magnitude of the peaks are found 
to be insensitive to temperatures kT <C A. 
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FIG. 3. The numerically calculated conductance as a func- 
tion of voltage for T L = 0.1, T R = 0.3, E Th = 0.4A, tf«A 
and N = 15. Inset: The differential conductance as a function 
of voltage for the diffusive SNS-junction, L = 0.29^tm and 
A = 0.17meV afc-fcT = 240mK, studied experimentally by 
Kutchinsky et alx3 The theoretical and experimental curves 
show the same qualitative features with conductance peaks at 
eV = 2A,A and A±E Th . 

These [features of the conductance have recently been 
observedli3£y in experiments with diffusive SNS-junctions 
of intermediate length, Exk % A (shown as inset in Fig. 

. The experimental curves show conductance peaks at 
eV w 2A,A,A±Exh and eV w (A±E Th )/2, with posi- 
tion and amplitude insensitive to temperatures kT <C A. 
In a corresponding interferometer setupEj the peaks at 
eV ks A ± Etk are demonstrated to result from induced 
PE (the peaks at eV ~ (A ± Exh)/2 could not be re- 
solved). The amplitude of all SGS peaks is successively 
reduced for increasing junction length, and the peaks 
at eV w A ± Exh are washed out in long junctions, 
Etk "C A. Although the system studied experimen- 
tally is different from a S-dot-S junction, the two type 
of junctions can be expected to show qualitatively simi- 
lar behavior, in the same way as for an NS and a N-dot-S 
junction (see Ref. |22| ). 

A qualitative understanding of the origin of the PE 
conductance peaks can be obtained by first considering 
the peak at eV w A + Exn ■ It is built up by quasiparti- 



cles which are Andreev reflected once at energies of order 
of Exh away from the Fermi energy (see Fig. |l|b). These 
scattering processes have enhanced amplitude due to the 
PE and are similar to the ones giving rise, to the finite 
bias conductance anomaly in NS-j unctions .□ The conduc- 
tance peak is insensitive to temperatures kT <C A, since 
the distribution of injected quasiparticles from the super- 
conductors is unaffected by temperatures well below the 
superconducting gap, in contrast to the NS-conductance 
anomaly, which is suppressed for kT 3> Exk ■ 

The conductance peak at eV « A — Etk results from 
scattering processes which include two Andreev reflec- 
tions, with one of the reflections at an energy of order 
Exh away from the Fermi energy. These process, in the 
same way, have an enhanced amplitudes due to the PE. 
Since the processes with one and two Andreev reflections 
are different, the shape and amplitude of the conduc- 
tance peaks at eV w A ± Exh are in general different. 
Following the same line of reasoning, we predict that 
the conductance (for Exh smaller than the distance be- 
tween subsequent subgap harmonics) will show peaks at 
alleF ^2{A ±E T h)/{2n). 

It follows from the quantum mechanical current ex- 
pression that the ensemble averaged dc-current can be 
written as (I) = N(e/h) + /£+ -/£" ), 

where fn h,± = {tr[{cj h,± y Cn /h ' ± }) fo/N (suppressing 
index a) are correlation functions of the wavefunction 
coefficients and the trace is taken over the transverse 
modes. The functions fn h,± can be interpreted as en- 
semble averaged distribution functions at energy E n for 
electron/hole quasiparticles with positive/negative sign 
of the wavenumber, summed over transverse modes. For 
broken time reversal symmetry in the dot, $ 3> <I> C , it is 
possible to formulate matching equations for the distri- 
bution functions fn^'^ directly, in the following way: 

We first note that for quasiparticles propagating in 
energy-position space along the MAR-ladder (see Fig. 
|l|b), the scattering matrices of the dot at different en- 
ergies are effectively uncorrelated to leading order in 
NTl,NTr, (i.e. neglecting quantum corrections). This 
result iS|-an extension of what is found for an N-dot-S 
junctionJlJ by employing the same diagrammatic tech- 
nique for integration over the unitary group. 

The statistical independence of the scattering matrices 
leads to three rules for wavefunction coefficient correla- 
tions and averages: i) coefficients with different energy 
indices n are uncorrelated since they are connected via 
at least one traversal through the dot. (This also means 
that there is no ac-Josephson current.) ii) coefficients 
for quasiparticles incoming towards the NIS-interface are 
uncorrelated, because incoming quasiparticles have scat- 
tered at dots at different energies before approaching the 
NIS-interface. iii) the average of any coefficient itself is 
zero. With these rules for the correlations between differ- 
ent wavefunction coefficients, we can, directly from the 
matching Eqs. (|1) and derive matching equations for 
the functions /„ ' . For the scattering across the dot, 
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e.g. for left injected electrons, we get from Eq. (||) 
f e n - = (tr [«'+ + t'c^) x c.c]) = l/2[#+ + f*'~ t ] 



f, 



n+1 



e,+l 



(5) 



In this derivation we used that the averaging rules gives 
Mc^W^-J) = (tr(^))(tr[(cS-)t^-])/JV and 
similarly for the other terms. Also, the averages 
(tr(t'tt')) = (tr(r+r)) = N/2. From Eq. (|), the match- 
ing equations at the left NIS-interface become 



te, 
J n 



f 

J n 



h,- 



= (tr [(b%(% + + + ctS n0 ) x c.c.]} 

= Bnfn + A n f n ' + + C n foS n o, 

= A n f*'+ + B n f%'- + D n f 5 n0 , (6) 



and 



where the scattering probabilities!!!]! A n 
similarly B ni C n and D n , have been introduced. The to- 
tal distribution function for right (left) going electrons on 
the left side at energy E' , f^,^_\(E') is found by sum- 
ming up the distribution functions for right(left) going 
electrons from all scattering states a. From Eq. (g), it 
follows that the total distribution functions at the left 
NIS-interface are related as 



ft{E')=A(E')[l-fl\ 
+ T(E')f (E'), 



-E')]+B(E')ft(E') 



(7) 



where T{E') = C(E') + D(E') = 1 - A(E') - B(E'). 
In this derivation we also used the symmetry for the to- 
tal distribution functions f e (E') = 1 — f h (—E'), which 
can be derived directly from Eqs. (H) and (|J). The re- 
lations between the total distribution functions at the 
right interface and across the dot follow in the same 
way. The resulting set of equations corresponds ex- 
actly to the matching equations for distribution func- 
tions in Ref. |l| (OTBK), with the important exception 
that the scattering by the dot itself couples the left- 
and right moving distributions of electrons in the nor- 
mal part of the junction (i.e. a SINFNIS junction with 
transparency 1/2 of 7'). The current is given by (I) = 
l/(eR )JdE'{f^(E>) - /«_(#)], with R = h/(2e 2 N). 
This shows that the OTBK-approach can be rigorously 
justified for a S-dot-S junction with broken time reversal 
symmetry in the dot. 

We notice that the presented derivation is also appli- 
cable to a long diffusive wire (/ <C L) SNS-junction with 
time reversal symmetry broken in the wire. The only 
difference is that the middle barrier (F) now has a trans- 
mission probability l/L. 

In conclusion, we have studied the dc-current transport 
in a voltage biased S-dot-S junction with an induced PE 
in the dot. It is found that the PE is manifested as peaks 
in the conductance at voltages eV « 2(A ± Exhl/^n. 
These peaks are insensitive to temperatures kT <C A but 
are suppressed by a weak magnetic field. The current for 
suppressed PE is independent of Etk and magnetic field 
and is shown to be given by the OTBK-theory. 



We thank J. Kutchinsky for making the experimental 
data presented in the paper available and we acknowledge 
helpful discussions with E. Bezuglyi, J. Kutchinsky, J. 
Bindslev-Hansen, J. Lantz and H. Schomerus. This work 
was supported by TFR, NEDO and NUTEK. 



1 M. Octavio et al, Phys. Rev. B. 27, 6739 (1983) 

2 A. Kastalsky et al., Phys. Rev. Lett. 67, 3026 (1991). 

3 A. A. Golubov and M. Yu. Kuprianov, Sov. Phys. JETP 
69, 805 (1989); W. Belzig, C. Bruder, and G. Schon, Phys. 
Rev. B. 54, 9443 (1996); S. Gueron et al., Phys. Rev. Lett. 
77 3025 (1996). 

4 J. Melsen et al., Europhys. Lett 35, 7 (1996). 

5 T.M. Klapwijk, G.E.Blonder, and M. Tinkham, Physica 
B+C 109-110, 1657 (1982). 

6 E.V. Bezuglyi et al., Phys. Rev. Lett. 83 2050 (1999), F. 
Pierre, ibid. 86, 1078 (2001). 

7 G.B. Arnold, J. Low Temp. Phys. 68, 1 (1987), E.N. Bra- 
tus', V.S. Shumeiko, and G.Wendin, Phys. Rev. Lett. 74, 
2110 (1995); D. Averin, and A. Bardas, ibid 75, 1831 
(1995); J.C. Cuevas, A. Martin-Rodero, and A.L. Yeyati, 
Phys. Rev. B. 54, 7366 (1996). 

8 N. van der Post et al., Phys. Rev. Lett. 73, 2611 (1994); 
E. Scheer et al., ibid. 78, 3535 (1997). 

9 A. Bardas and D. Averin, Phys. Rev. B. 56, R8518 (1997). 

10 Y. Naveh et al., Phys. Rev. Lett. 85, 5404 (2000). 

11 E.V. Bezuglyi et al., Phys. Rev. B 62, 14439 (2000). 

12 B.A. Aminov, A. A. Golubov, and MYu. Kupriyanov, 
Phys. Rev. B. 53, 365 (1996), A.V. Zaitsev and D.V. 
Averin, Phys. Rev. Lett. 80, 3602 (1998), E. Scheer et al. 
ibid 86, 284 (2001). 

13 J. Kutchinsky et al., Phys. Rev. B 56, R2932 (1997); R. 
Taboryski et al., Superlatt. and Microstr. 25, 829, (1999). 

14 See e.g. G. Johansson et al. Superlatt. and Microstr. 25, 
906, (1999) and references therein. 

15 C.W.J. Beenakker, Rev. Mod. Phys, 69, 731 (1997); Y. 
Alhassid, ibid 72, 895 (2000). 

16 G.E. Blonder, M. Tinkham and T. M Klapwijk, Phys. Rev. 
B. 25, 4515 (1982). 

17 P.W. Brouwer and C.W.J. Beenakker, J. Math. Phys. 
(NY.) 37, 4904 (1996). 

18 The excess current is ranging from (Gjve/A)[16/3+7r(9/2— 
Ay/2)] for T = 1 to (Give/A)[7r/4(7 - 4^2)] for T < 1. In 
the limit r <JC 1 the distribution p(D) is the same as for a 
dirty tunnel barrier and the result coincides with Ref. [To| 

19 In all numerical curves shown, M — 20N and 1000 Hamil- 
tonians have been generated for each curve. 

20 For large interface transparency rz,,r.R ~ 1, the conduc- 
tance peaks turn into dips. 

21 Additional features in the SGS has also been observed in A. 
Chrestin, T Matsuyama and U. Merkt, Phys. Rev. B 55, 
8457 (1997) and T. Schapers et al, Superlatt. and Microstr. 
25, 851, (1999). 

22 A. A. Clerk, PW. Brouwer, and V. Ambegaokar, Phys. Rev. 
B 62, 10226 (2000). 



4 



